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Abstract. We present recent lattice QCD results on the electroweak nucleon to A transition and A 
form factors using dynamical fermion gauge configurations with a lowest pion mass of about 300 
MeV, with special emphasis in the determination of the sub-dominant quadrupole Nf — > A and A 
electromagnetic form factors. 
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INTRODUCTION 

There is a number of recent lattice QCD calculations of the electroweak form factors of 
the nucleon [1]. The main focus of this presentation is the evaluation of the electroweak 
nucleon (N) to A and A form factors (FFs). In Ny* — > A the dominant magnetic dipole FF, 
G* m , is precisely measured and it therefore serves, like in the case of the electromagnetic 
(EM) FFs of the nucleon, as a benchmark of the lattice QCD methodology. The sub- 
dominant quadrupole FFs G* E2 and G* C2 have also been studied extensively, since their 
value carries information on the deformation in the N/A system. In order to calculate 
them in lattice QCD, one applies several improvements to attain good enough accuracy. 
The dominant axial N to A form factors Cg and Cg can also be calculated within lattice 
QCD. They correspond to the nucleon axial, Ga, and nucleon induced pseudo-scalar, 
G p , FFs and provide important input for phenomenological models and chiral effective 
theories. Furthermore, the evaluation of the kNA coupling enables us to check the 
validity of PCAC and the associated non-diagonal Goldberger-Treiman relation. 

The FFs of the A are difficult to measure experimentally due to the A short lifetime. 
Only the magnetic moment of the A is measured albeit with large uncertainty. Thus 
lattice QCD can provide a valuable input on the A FFs enabling for instance the deter- 
mination of the A charge distribution in the infinite momentum frame [2]. Due to the 
3/2-spin structure of the A, there are two Goldberger-Treiman relations, which can be 
examined by calculating, besides the axial A FFs the pseudo-scalar ones. Having the 
EM and axial FFs for the N/A system enables a combined chiral fit to determine the 
coupling constants that enter in chiral effective models. 

Recently there has been a lot of progress in the calculation of hadron masses and 
nucleon structure using dynamical lattice QCD simulations [3]. The computational cost 
of these simulations can be parameterized as a function of the lattice spacing a, the 

lattice spatial extent L and the pion mass, «%, as C sim oc ( 30 ™ eV ) ™ (jffi)^ (^rY" ■ 
The coefficients c m , cl and c a depend on the type of discretized action. Based on current 



simulations the cost at the physical point is estimated of &{\) Teraflop-year. 
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FIGURE 1. Left: The low-lying baryon spectrum computed using Nf = 2+1 Clover [4] and Nf = 2 
twisted mass fermions (TMF) [5]; Upper right: The nucleon axial charge [6]; Lower right: The nucleon 
isovector moment of the unpolarized quark distribution [7]. The physical point is shown by the asterisk. 

In Fig. 1 we show recent results on the low-lying baryon spectrum obtained by the 
BMW Collaboration using Nf = 2 + 1 Clover fermions and the ETM Collaboration 
using Nf = 2 twisted mass fermions. Both collaborations used 3 lattice spacings to ex- 
trapolate the results to the continuum limit. One observes that the results using different 
discretization schemes are in agreement and that both reproduce the experimental val- 
ues. This is a significant validation of lattice QCD techniques. 

The electric G£-(g 2 ) and magnetic Gm(<] 2 ) Sachs form factors (FFs) as well as the 
axial- vector FFs Ga(<? 2 ) and G p (q 2 ) of the nucleon have also been studied by sev- 
eral lattice groups using dynamical simulations down to lowest pion mass of typically 
~ 250 MeV. Lattice data are in general agreement, but still show discrepancies with ex- 
periment. In Fig. 1 we compare recent results from various collaborations on the nucleon 
axial charge and isovector moment of the unpolarized quark distribution. Whereas there 
is an overall agreement among lattice results the experimental values are not reproduced 
(see Ref. [1] for more details). 



Nf ->• A TRANSITION FORM FACTORS 

The Ny* — > A transition is written in terms of three Sachs FFs: 



(A(p',s')\j^\N(p,s)}= £ /u a {p',s' 



with srf = zy§ ( e a (&)En(p) ) • There is a wealth of experimental information on 
the Ny* — > A transition [8]: The dominant magnetic dipole FF, G* Ml , is well mea- 



G *Ml (l 2 ) K ail + G *E2(l 2 ) K al + G C2 K al 



u(p,s) (1) 



sured and the electric G E2 and Coulomb G£ 2 FFs are found to be non-zero signal- 
ing a deformation in the nucleon/A-system. The deformation is probed via the ratios: 

Rem (EMR) = - §Pjp , and R SM (CMR) = - i tIHtSt > in tne rest frame of the A - As 



shown in Fig. 2, precise data strongly "suggest" deformation of the N and/or A [9, 10]. 
New data on CMR at low momentum transfer are currently being analyzed [11]. 
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FIGURE 2. Left: The ff Lr of the p(e,e'p)n° reaction at Q 2 = 0.127 GeV 2 [9]; Right: Experimental and 
lattice QCD results on the magnetic dipole FF G* Ml as a function of Q 2 = —{p' — p) 2 - 

The lattice evaluation involves the computation of two- 
point and three-point functions: 




(*„«,) 



G(q,t) = £e-"/'?r 4 (J h (x f ,t f )J h (0)) 
G» v (r,q,t) = J^j^r(J h (x f ,t f )^(x,t)J h (0)) 

Xf,X 

where Jh{x,t) is an interpolating field creating a state with the quantum numbers of the 
baryon h and we have taken tj = 0. One computes the three-point function at various 
^-values of the current insertion which, for large tt and t and taking an appropriate 
ratio with two-point functions, yields the matrix element of Eq. (1). 
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FIGURE 3. Results on EMR (left) and CMR (right). 

In order to extract the sub-dominant quadrupole FFs one constructs optimized sources 
to isolate them from the dominant dipole and uses the coherent sink technique to 
increase statistics. Recent lattice results are shown in Figs. 2 and 3 using a hybrid action 
of dynamical staggered sea and domain wall valence quarks, as well as Nf = 2 + 1 



dynamical domain wall fermions (DWF), simulated by the RBC-UKQCD with lowest 
pion mass of about 300 MeV [12]. The slope of G* Ml at low Q 2 remains smaller than 
what is observed in experiment underestimating G* Ml (0) i.e. one observes the same 
effect as for the nucleon form factors. Since G* E1 and G* C2 are underestimated at low 
Q 2 like G* Ml taking ratios may remove some of these discrepancies. Indeed the EMR 
shown in Fig. 3 is in better agreement with experiment, whereas CMR approaches the 
experimental values as the pion mass is lowered. Despite the increased statistics the 
errors on the sub-dominant ratios are large and to reduce the errors as m K approaches 
its physical value one would need to increase significantly the number of statistically 
independent evaluations. 



N - A AXIAL-VECTOR AND PSEUDO-SCALAR FORM FACTORS 

The N - A axial- vector matrix element (&{p' ,J)\AV\N{p,s)) is written as 



m N 



f + 



P ,V {gXiigpv ~ glpgfiv) q P ' +C^{q 2 )g 



-qxqfi 



whereas the N - A matrix element of the pseudo-scalar current is given by 



2m q (A(p' lS ')\P 3 \N(p lS )) 



J K m 2 GjtNAjq 2 
m l % -q l 



-u v (p ,s)- u(p,s 

2mjy 



u{p,s), (2) 



(3) 



Using the axial Ward identity and pion pole dominance one obtains the non-diagonal 
Goldberger-Treiman (GT) relation, G^^q 2 ) f n = 2m^C^(q 2 ). In Fig. 4 we show results 
on the dominant axial FF Cg and on the ratio G^A(2 2 )A/2mA?Cg (Q 2 ), which should 
be unity if the GT relation holds. This ratio approaches unity for Q 2 > 0.5 GeV 2 [12]. 
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FIGURE 4. The dominant axial and pseudo-scalar N to A FFs using a hybrid action and DWF. 



A FORM FACTORS 



The A matrix element of the electromagnetic current (A(p ; ,s') \j^ \A(p, s)) is given by 



with e.g. the quadrupole FF given by: Gei = (F* — tF 2 * ) — i(l + t) (F 3 * — tF 4 * ), where % = 
—q 2 /(4Ml). Using lattice results on Gei one can obtain the transverse charge density of 
a A in the infinite momentum frame [13, 14]. This is shown in Fig. 5, where a A with 
spin 3/2 projection along the x-axis is elongated along the spin axis [13]. In the same 
figure we also show the corresponding charge density of the Q. , which shows a similar 
deformation as the A [15]. 
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FIGURE 5. Left: Lattice results on the A electric quadrupole FF; Contours of A (middle) and Q. 
(right) [15], with 3/2 spin projection along the x-axis. Dark colors denote small values. 

The A matrix elements of the axial-vector current (A(//, s') \A\ |A(p, s) ) is given by 

-1 



Ua(p ,S 



up{p,s) (5) 



and of the pseudo-scalar current (A(p f ,s')\P i \A(p,s)} by 



-u a {p',s' 



2m q (ml-q 2 ) 



up(p,s) 



(6) 



The A axial charge is derived from gi(0) [16, 17], whereas there are two 7zAA pseudo- 
scalar FFs. 
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FIGURE 6. Lattice QCD results on the dominant axial and pseudo-scalar A FFs in the quenched theory 
and using a hybrid action [17]. 

G^aa(0) is non-zero and can be identified as the % — A coupling. Lattice QCD results 
on the dominant axial A FFs g\ and #3 as well as on G^aa are shown in Fig. 6. One can 
derive two Goldberger-Treiman relations: fnG^^q 2 ) = m Agi(l 2 ), and fnH^^q 2 ) = 
m^h\ (q 2 ), which can be tested using lattice QCD results. 



CONCLUSIONS 



We have shown that lattice QCD successfully reproduces the low-lying baryon spectrum 
using different discretization schemes. There is an on-going investigation of nucleon 
structure by a number of lattice collaborations. Similar techniques can be applied to 
study transitions and resonant properties and we have applied these methods in the study 
of the N to A electroweak transition FFs as well as the A FFs. The latter are difficult 
to measure experimentally and therefore lattice QCD provides valuable input on these 
quantities. Having lattice QCD results on the N-A system one can use, for the first time, 
chiral perturbation theory to extract the axial couplings gA, ca and g& from a combined 
chiral fit to the lattice results on the nucleon and A axial charges and the axial N to A form 
factor Cs(0) [17]. Applying such a fit to lattice results in the pion mass range from 500 
MeV to 300 MeV, still does not reproduce the experimental value of gA- Current lattice 
QCD simulations reaching pion masses below 200 MeV are now becoming available 
and these simulations, combined with a detailed study of lattice systematics [18], are 
expected to shed light on the origins of the observed discrepancies. 
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